Theory of microwave-assisted supercurrent in quantum point contacts 
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We present a microscopic theory of the effect of a microwave field on the supercurrent through a 
quantum point contact of arbitrary transmission. Our theory predicts that: (i) for low temperatures 
and weak fields, the supercurrent is suppressed at certain values of the superconducting phase, (ii) 
at strong fields, the current-phase relation is strongly modified and the current can even reverse 
its sign, and (iii) at finite temperatures, the microwave field can enhance the critical current of the 
junction. Apart from their fundamental interest, our findings are also important for the description 
of experiments that aim at the manipulation of the quantum state of atomic point contacts. 

PACS numbers: 74.40.Gh, 74.50.+r, 74.78.Na 



The DC Joscphson effect is one of the most striking 
examples of macroscopic quantum coherence. In the con- 
text of superconductivity this phenomenon manifests as 
the flow of a dissipationless DC current through a junc- 
tion in the absence of any voltage [l], Q • Since its discov- 
ery in the early 1960s, this effect has been observed in a 
large variety of weak links such as tunnel junctions 
microbridges atomic contacts carbon nanotubes 
0, semiconductor nanowires 0] and graphene Q. In 
spite of their intrinsic differences, the DC Josephson ef- 
fect in these systems can be described in a unified man- 
ner. It has been shown that for constrictions shorter 
than the superconducting coherent length, the Joseph- 
son current is carried exclusively by a single pair of An- 
dreev bound states (ABSs) In the simplest case of 

a single-channel contact of transmission r, these states 
appear at energies E^{ip,T) — ±A[1 — r sin^((y5/2)]-^/^, 
where A is the superconducting gap and (p is the phase 
~^ ' difference between the order parameters on both sides 
I (see Fig. [T]) . These states carry opposite supercurrents 
^ . /±(<p) = {2e/h)dE^/dip, which are weighted by the oc- 



spectroscopy. This leads us to the central question of 
the present work: How does a microwave radiation mod- 
ify the supercurrent of a single-channel quantum point 
contact (QPC)? Apart from its fundamental interest for 
the field of mesoscopic superconductivity, this question 
is also of great relevance for the proposals of using the 
ABSs of a QPC as the two states of a quantum bit [12I - 
14l |. whose quantum state can be probed by means of 
current measurements. 

Surprisingly, there is no complete answer to the ques- 
tion posed above. The theoretical analysis of the 
microwave-assisted supercurrent in point contacts has ei- 
ther been addressed within phenomenological approxi- 
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cupation of the ABSs. A more complex weak link, like 
the ones mentioned above, can be viewed as a collec- 
tion of independent conduction channels, characterized 
by a set of transmission coefficients. The supercurrent 
through it is given by the sum of the contributions from 
the individual channels [IQ] . 

This unified microscopic picture of the DC Josephson 
effect has been recently confirmed experimentally in the 
context of atomic contacts , where the current-phase 
relationship has been directly measured. These experi- 
ments mainly probed the ground Andreev state and, as it 
is stated in Ref. [nj , it would be of great interest to probe 
also the excited state, for instance, through microwave 
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FIG. 1. Dispersion relation of the ABSs in a single-channel 
QPC with transmission r = 0.95. The vertical dashed lines 
indicate one- and two-photon transitions between the ABSs 
(a-b), and transitions between the continuum and ABSs (c-d). 
We have chosen i? = at the Fermi energy. 
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mations or in the limiting case of very weak fields 
In this letter, we present a microscopic theory of the ef- 
fect of a microwave field on the supercurrent of a single- 
channel quantum point contact valid for arbitrary range 
of parameters. Our theory based on the Keldysh tech- 
nique predicts the following novel effects: (i) at low tem- 
peratures, the supercurrent can be strongly suppressed at 
certain values of the phase due to resonant microwave- 
induced transitions between the two ABSs (processes of 
type a and b in Fig. [1]). (ii) As the radiation power in- 
creases, the supercurrent-phase relation is strongly modi- 
fied and it can even reverse its sign, (iii) At finite temper- 
atures, the radiation can induce the transition of quasi- 
particles from the continuum to the lower ABS leading 
to an enhancement of the critical current as compared to 
the case in the absence of microwaves (process of type 
c in Fig [T]) . We also compare our results to a two- level 
model (TLM), where the QPC is described exclusively in 
terms of the ABSs and show that effects (ii) and (iii) fall 
out of the scope of TLM. This is especially relevant for 
the quantum computing applications. 

We consider a QPC consisting of two identical super- 
conducting electrodes (denoted as L and R) and linked 
by a single conduction channel of transmission r. Our 
goal is to compute the supercurrent through this QPC 
when it is subjected to a monochromatic microwave field 
of frequency w. We assume that the external radiation 
generates a time-dependent voltage V{t) = Vbsinw^ 0. 
According to the Josephson relation, this voltage induces 
a time-dependent superconducting phase difference given 
by (/)(i) = (fi + 2a cos w^, where ip is the DC part of the 
phase and a = eVo/huj is a parameter that measures the 
strength of the coupling to the electromagnetic field and 
that is proportional to the square root of the radiation 
power at the junction. Following Refs. [l7, 18|, the cur- 
rent through a QPC with an arbitrary time-dependent 
voltage can be computed as I{t) = (e/4?i)Trf3/^(i, t), 
where T3 is the third Pauli matrix and I^{t,t) is the 
Keldysh component of the current matrix given by 

i{t, t')=2T[GL,Gj,]^o[4~T{2~{GL,GR}jy\t,t'). 

(1) 

Here the symbol " represents 4x4 matrices in Keldysh- 
Nambu space and the symbol o denotes the con- 
volution over intermediate time arguments. More- 
over, Gl{r) are the quasiclassical Green functions 
for the left and right electrodes, which can be ex- 
pressed as Gj{t,t') = e''^^(*)^3/2^^.(i _ i')e-''^^(*')^^/2^ 
Here, g{t) = J {dE /2TT)e-''^*/^g{E) is the equilibrium 
Green function of the leads and (t) is the time- 
dependent phase of the j superconductor, j = L,R, 
i.e. 0L(i) ~(t)Bit) = 0(i)/2. The retarded (i?), 
advanced {A) and Keldysh {K) components of g{E) 
adopt the form: g^^^^^E) = g^^^'^{E)Ty. + ]^^^'^{E)iT2 
and g^{E) = {g^[E) - g^{E)\i&^\i{E jlk^T), where 
gRi^){E) = E/[{E±i7^f - l^^y/^, f^'^^\E) =/^/[{E± 
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It is easy to show that, due to the time dependence 
of the phase, the lead Green functions Gl{r), and any 
product of them, admit the following Fourier expan- 
sion Git,t') = E::=-ooe"""*7ife"'''^*"*'^^''Go™(i?), 
where Gnm{E) = G{E + nhuj, E + mhuj) are the corre- 
sponding Fourier components in energy space. Thus, / in 
Eq. ([1]) can be written as a product of matrices in energy 
space. In particular, its Keldysh component is given by 
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Here, we have defined the matrices Anm = 2t[Gl, Gii\nm 
and = [41 - t(2 - {Gl, Gi^})]"^, which can be 

determined from the Fourier components of G]^(^fjy Once 
the components of are obtained from Eq. ^ , one can 
compute the current. We are only interested in the DC 
component, which reads 
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Trf3/o^(£;,(p,^,a). (3) 



The DC current can only be calculated analytically in 
certain limiting cases like in the absence of microwaves, 
in the tunnel regime or for very weak fields. In general, 
Eq. ^ and the current have to be evaluated numerically. 

In the absence of microwaves, the current from Eq. (|3]) 
can be written as a sum of the contributions of the two 
ABSs as /oq('^) = IJ^fiFiEJ^) + I^n-FiEX), npiE) being 
the Fermi distribution function, which yields 



sA^ rsiniy9 
'2hE+{^ 



tanh 



2kBT 



(4) 



In equilibrium, the maximum current is then obtained 
at temperature T = 0, when the lower ABS is fully oc- 
cupied and the upper one is empty. In the presence of 
the microwave field the siinplest approach is the so-called 
adiabatic approximation [3|, which consists in replacing 
the stationary phase ip in Eq. (j4]) by the time-dependent 
phase (pit). This leads to 



hdi'P, a) = X! InJo{2na) sin(n<^). 



(5) 
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where /„ — (I/tt) J^^ dip Icq{p) sm{np) are the harmon- 
ics of the equilibrium current-phase relation and Jq is the 
zero-order Bessel function of the first kind. 

In Fig. [2][a) we show the zero-temperature current- 
phase relation (CPR) computed numerically from Eq. (|3]) 
for a highly transmissive channel with t ~ 0.95 and a 
weak field a = 0.1 of frequency hio = 0.6A. For com- 
parison we also show as a dashed line the result obtained 
with the adiabatic approximation of Eq. ([5]). The main 
difference is that the exact result shows a series of dips 
where the current is largely suppressed. These dips orig- 
inate from microwave-induced transitions from the lower 
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FIG. 2. (Color online) Zero-temperature supercurrent, in 
units of lo = eAo/h, as a function of the phase, (a) The 
solid line shows the exact result for r = 0.95, ho = 0.6A and 
a = 0.1. The dashed line corresponds to the approximation 
of Eq. ([5]) and the dotted line to the TLM. Rest of the panels 
show the exact result for (b) r — 0.95, a = 0.1, and different 
frequencies, (c) r = 0.95, hoj = 0.3A, and different a's, (d) 
huj = 0.3A, a = 0.1 and different r's. 



ABS to the upper one that enhance the population of the 
latter, diminishing the supercurrent (see processes a and 
b in Fig. [1}. Such transitions can occur whenever the 
Andreev gap (distance between the ABSs) is equal to a 
multiple of the microwave frequency, i.e. 2E^{ip) — nhxjj, 
where n = 1,2,... can be interpreted as the number of 
photons involved in the transition. For small values of 
a the ABSs remain almost unchanged, thus the resonant 
processes take place at phases given by 



2arcsinv/[l - {nhLj/2A)'^]/T, n = 1, 2, 



(6) 



This expression reproduces accurately the positions of 
the dips in Fig.[2ja). The origin of the dips can be further 
confirmed by exploring the CPR for different frequencies, 
as we do in Fig. [51Jb). Here, one can see that by decreas- 
ing the frequency, the dip of order n = 1 moves to higher 
values of (p and disappears for fku < 0.4 A, in agreement 
with Eq. dH). 

One can gain further insight by analyzing this problem 
in terms of a TLM that describes the dynamics of our 
QPC in terms of the ABSs [H, [H. We consider the 



TLM of Ref. [13[, whose effective Hamiltonian in the 
instantaneous basis of ABSs reads 



HA{t) ^ E+{^{t))a, 



rrA2 sin2(0(i)/2) 



mt)cTy, (7) 



where ay^z are Pauli matrices, r = \/l — t and 0(i) = 
d(t>{t)/dt. We have computed the CPR from this model 
using a Floquet approach 20]. In Fig. we show a 
comparison of the results of this TLM with the exact 
results. There is an excellent agreement in a wide range 
of phases and in particular, the TLM is able to reproduce 



the current dips. A discrepancy occurs at phases close 
to TT, which is understandable as the model assumes that 
h(t){t) ^ atu ^ 2£'^, while for (p ^ tt and high r, the 
ABSs are very close to each other and the assumption no 
longer holds. Using a rotating wave type approximation 
we also obtain from the TLM the following analytical 
expression for the DC current at the first two resonances. 
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Here and below. 



and E'l 



are the first and second 



derivatives of E^ with respect to the phase. The detun- 
ings Si 2 are given by i5i = E^{ip) — huj/2 + a'^E'l{ip) + 
3r2a2fil^(A2-£;+(<^)2)2/(32£;+((^)4), ^nd 82 = E+{ip) - 
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include shifts in the ABS energies caused by the mi- 
crowave field. The resonance widths are 71 — rafta;(A^ — 
El{ipY)/{4.E+{>pf) and 72 = ra^hojEl{ip){2A^ - 
E\{ipY)/{2E^{lpY), in general being proportional to the 
coupling strength, 7„ oc a". For small power, Eq. ([SJ is in 
a good agreement with the exact solution for the current- 
phase relation. We also find that within the TLM the 
current vanishes completely at the resonances (given by 
the condition (5„ = 0), as a consequence of the fact that 
at these points the long-time average populations on the 
ABSs are equal, so that the currents carried by the two 
states cancel exactly 2l|. Transitions from the contin- 



uum, however, make the current at the resonances finite 
and dependent on the frequency and the excitation ener- 
gies. This is seen for example in the exact result for the 
second dip shown in Fig. [5] ^||. 

Let us now discuss the dependence of the CPR on 
the radiation power. In Fig. [^c) we show the CPR for 
T = 0.95, hhj — 0.3A and different values of a. As a in- 
creases, the CPR is drastically modified and the current 
is not only strongly suppressed around the phase values 
given by Eq. but everywhere. Notice also that in 
certain regions, particularly at large phases, the current 
even reverses its sign. It is also worth stressing that as 
a increases, we find larger deviations between the ex- 
act results and those of the TLM (not shown here) due 
to multi-photon processes connecting the ABSs and the 
continuum of states. For the sake of completeness, we 
illustrate in Fig. [DJd) the influence of the transmission in 
the CPR for hw — 0.3A and a — 0.1. In this case, for 
T < 0.9 the CPR can be accurately described with the 
adiabatic approximation for all phases since the transi- 
tions between the ABSs are very unlikely. 

We turn now to the analysis of the critical current Ic^ 
i.e. the maximum value of the DC Josephson current. In 
Fig.IS^a) we show the critical current as a function of a for 
several temperatures. One clearly sees that the adiabatic 
approximation (dashed lines) only describes correctly the 
behavior of Ic when the frequency, power and tempera- 
ture are low enough so that the microwaves cannot induce 
transitions between the ABSs and between them and the 
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FIG. 3. (Color online) Panels (a-b) show the critical current as a function of a for (a) hw = 0.3A, r = 0.95 and different values 
of T, and (b) f\u) — 0.3A, ksT — 0.5A and different values of r. In both panels the solid lines correspond to exact results and 
the dashed ones to Eq. ((5]). (c) Current as a function of the frequency for ip = 2.0, kuT = 0.5A, t = 0.95 and a = 0.1. 



continuum. The most striking result in Fig. ^a) is the 
enhancement of Ic with respect to the case a = as 
the temperature is raised. In Fig.[3][b) we show that this 
occurs in a wide range of transmission values. The en- 
hancement of the critical current by a microwave field has 
been predicted and observed in superconductin g in icro- 
bridges 22 , 2^ and in proximity effect structures 13, 2^ ■ 
In both cases the enhancement is due to a redistribution 
of the excitations induced by the field. In our case the 
underlying mechanism is similar. At finite temperature 
the lower ABS is not fully occupied and the microwave 
field can promote quasiparticles from the continuum to 
this state if huj > A + E^{(p), for values of ip close to that 
where the maximum takes place (see process c in Fig.[T]). 
This naturally results in an enhancement of the super- 
current. We illustrate this argument in Fig. ISjc), where 
we show the frequency dependence of the supercurrent 
for a fixed value of the phase (po = 2.0. As one can see, 
the current remains largely unaffected until hu reaches 
A -I- E^{(po) Ki 0.46A, where it starts to increase. Then, 
as the frequency increases further one can observe the ap- 
pearance of two dips, corresponding to the resonances at 
Huj = 1 2 and E\ . For larger values of ftw one finally 
sees a decrease of the current at A + E'^iip^) « 1.54A 
due to processes that promote quasiparticles from the 
continuum to the upper ABS (see process d in Fig. [ij. 
Due to electron-hole symmetry similar transitions exist 
between the ABS and the upper continuum. Obviously, 
the enhanced supercurrent cannot be explained in terms 
of any TLM since it involves the continuum of states. 

In summary, we present here a microscopic theory of 
the microwave-assisted supercurrent in quantum point 
contacts. It predicts the appearance of a variety of novel 
phenomena that in general are out of the scope of simple 
approximations and two-level models. Our results are 
of relevance for many different types of weak links and 
in particular, they can be quantitatively tested in the 
context of atomic contacts. 
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